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Abstract 

We calculate the 0(a s ) and 0(af) gluon radiative corrections to the QCD sum 
rule for the first Gegenbauer moment af of the kaon light-cone distribution ampli- 
tude. The NNLO accuracy is achieved for the perturbative term and quark-condensate 
contributions to the sum rule. A complete factorization is implemented, removing 
logarithms of s-quark mass from the coefficients in the operator-product expansion. 
The sum rule with radiative corrections yields a^(l GeV) = 0.10 ± 0.04. 
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1. Light-cone distribution amplitudes (DA's) of hadrons enter various factorization 
formulae used for description of exclusive processes in QCD. The concept of DA's allows to 
describe collinear partons in an energetic hadron, separating long-distance dynamics from 
the perturbatively calculable hard-scattering amplitudes. 

The set of DA's with a growing twist is especially useful for the pion and kaon be- 
cause their intrinsically small masses make the collinear description more efficient. The 
lowest twist-2 DA has a transparent physical interpretation, describing the longitudinal 
momentum distribution in the quark- ant iquark Fock-state of a meson. Switching from the 
pion to kaon, one encounters the SU(3) /^-symmetry violation effects, which originate from 
the quark mass difference m s — m u ^- These effects have to be accounted as accurate as 
possible, in order to assess the SU(3)fi symmetry relations between the hadronic ampli- 
tudes with pions and kaons. Important examples are the relations between B — > -kit and 
B — > ttK, KK charmless decay amplitudes employed in the studies of CP- violation and 
quark-flavour mixing. 

The most essential SU(3) fi -violating effects in the kaon twist-2 DA include the ratio of 
the decay constants Jk/ fn and the difference between the longitudinal momenta of strange 
and nonstrange quark-partons. This difference is proportional to the first moment of in 
the decomposition of the kaon twist-2 DA in Gegenbauer polynomials, whereas vanishes 
in the isospin (G-parity) symmetry limit. In addition, the ratio of the second Gegenbauer 
moments af /cig can also deviate from unity; the effects related to a% at n > 3 are usually 
neglected. 

In this paper we concentrate on the determination of the asymmetry parameter af (/i) 
for the kaon, at a low scale \x ~ 1 GeV. The method originally suggested in [1] and 
based on QCD sum rules p] is employed. The most recent sum rule estimates of af were 
obtained in [3] and pE], where, in addition to the known leading-order (LO) results, the 
next-to-leading (NLO), 0(a s ) correction to the quark-condensate contribution are taken 
into account. These calculations, together with the estimates [5], [6] based on the operator 
identities, yield the interval (quoted as a best estimate in [7]): af (lGeV) = 0.06 ± 0.03. 
The positive sign of af corresponds, as expected, to a larger average momentum of the 
heavier valence s-quark in the kaon. 

The aim of this work is to upgrade the precision of the QCD sum rule for af . We 
calculate the gluon radiative corrections to the perturbative and quark-condensate contri- 
butions in NNLO, including 0(a s ) and 0(a 2 ) corrections. This task is technically feasible, 
due to the currently achieved state-of-the-art in the calculations of multiloop effects in the 
two-point correlation functions with strange and nonstrange quarks. For the correlation 
function with the scalar and pseudoscalar currents the O(af), five-loop accuracy has re- 
cently been achieved [8], and used, e.g., for the QCD sum rule determination of the strange 
quark mass 0, [M [EE]- In this case, in the perturbative expansion the 0(a 2 ) terms are 
important numerically, which is one motivation to include these terms also in the sum rule 
for af" . The correlation functions underlying the sum rules for Gegenbauer coefficients are 
however different, because the currents contain derivatives. Therefore, the calculation re- 
ported in the present paper, involves a certain technical novelty. In addition, we clarify and 
take into account the mixing of operators that is necessary for the complete factorization 
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of small and large scales in the correlation function. Our result for the first Gegenbauer 
moment of the kaon DA is: 

of (1 GeV) = 0.10 ±0.04. (1) 

In what follows, we introduce the correlation function, present the expressions for the 
new radiative corrections, derive the resulting QCD sum rule for af , including the new 
0(a s ) and 0(a 2 ) terms and perform the numerical analysis. 

2. The twist-2 DA of the kaon enters the standard expression for the light-cone expan- 
sion of the vacuum- kaon bilocal matrix element (we take K~ for definiteness) : 

(K-(q)\s(z) llMl5 [z, -z] u(-z)\0) z ^ = -tq,f K f du e iu ^ i ^ z < fK (u, fi) , (2) 

Jo 

where the s- and u quarks carry the momentum fractions u and u = 1 — u; [z, —z] is 
the path-ordered gauge-factor [xi,X2] = Pexp(i J* dv(x\ — X2) p A p {vxi + VX2)) , and fi is 
the normalization scale determined by the interval z 2 near the light-cone. We use the 
compact notation A p = g s A a p \ a /2 for the gluon field and the covariant derivative is defined 
as Dp = d p — iA p . In (T5]), the twist-2 DA ipx{u) is normalized to unity, so that in the local 
limit z — > one reproduces the definition of the kaon decay constant fx- 
As usual, (Pk(u) is expanded in the Gegenbauer polynomials 

y K {u,n) = 6uu (l + |>*(//)C'f 2 ( U - u)J , (3) 

with the coefficients a^(fi) (Gegenbauer moments). The first Gegenbauer moment of is 
proportional to the average difference between the longitudinal momenta of the strange and 
nonstrange quarks in the two-parton state of the kaon. Expanding both parts of Eq. (j2j) 
around z = in local operators and using the decomposition (j3J) with C^ 2 (x) = 3x one 
relates af to the vacuum-to-kaon matrix element of a local operator with one derivative: 

<-» 3 
(K~{q)\s r Yv r y5iD\u\0) = -iq v q x f K -= a x 1 ( 4 ) 

5 

«— > — > * — 

where D\=D\ - D\. 

The Gegenbauer moments a^ ,K (fi) are known to be multiplicatively renormalizable 
only at the one-loop level. Generally, this property is lost at higher orders in a s , e.g., 
the two-loop renormalization of calculated in [12] includes operator-mixing effects. Still 

the af case is special, in so far as the underlying operator 57^752 D\ u can only mix with 
d\{s^ v 75«), as there is no other local operator with the same dimension and flavour content. 
However, the above two operators have opposite G( s )-parities, where G/ 8 ) is the analog of 
the isospin G-parity for the 577(2) subgroups of SU(3)fi involving s quark {V or [/-spins). 
Naturally, G( s )-conservation is only realized in the m s = m u ^ limit. Note, however, that 
the ultraviolet renormalization in MS-scheme is a mass-independent procedure. Hence it 
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is legitimate to consider the SU(3)fi limit, while performing the renormalization, so that 
the G( s )-conservation protects the operators from mixing with each other. As a result, af 
remains multiplicatively renormalizable at any order in perturbation theory. For complete- 
ness, we present the well-known expression for scale-dependence of af with the two-loop 
(NLO) accuracy, written in an unexpanded form: 

7Q / 7Q 71 \ 

a*W =(^m" (§±fp<^M U " SJ a*M , (5) 
\a a {fio)J V Po + Px[a 8 {n)/it) J 

where 7 = 8/9, 71 = 590/243 are the anomalous dimensions [T3] and (3q — 9/4, /5i = 4 
are the coefficients of /3-function for rif — 3. 

As originally suggested in [1], the few first Gegenbauer moments of DA's can be cal- 
culated employing operator-product expansion (OPE) and QCD sum rules for two-point 
vacuum correlation functions. The method works well only for the first two coefficients 
al'2 ■ In the sum rules for a^y 3 the condensate contributions grow fast with n and the 
control over OPE is lost H 

To obtain a QCD sum rule for af, it is convenient to take the so called "diagonal" 
correlation function: 

IW?) = *jd 4 x e^'(0|T ^u(x)^ l5 s(x), s(0) 7 „75* Sa«(0)} |0> = q^qx^q 2 ) + ■■■, (6) 

where, for brevity, only the relevant kinematical structure is shown. In fact, (EJ) is not quite 
diagonal, because one of the operators, the same as in (J3J), contains a derivative. A differ- 
ent choice is to correlate the operator in (j3J) with the pseudoscalar current ui^^s. This, so 
called "nondiagonal" correlation function was tried in [15], but produces an unstable sum 
rule (see the discussion in [5]). 

3. In LO, the operator product expansion (OPE) for the correlation function ([6]) 
includes the quark-loop diagram at 0(a°) and the contributions of the vacuum condensates, 
calculated in the deep spacelike region, Q 2 = —q 2 ^> Aq CD . The corresponding diagrams 
are collected in Fig. [TJ One obtains a generic expansion for the invariant amplitude II(g 2 ) 
defined in OH]) in inverse powers of the variable Q 2 : 

wr) 2 , x _ MQ\ri , A(Q 2 ,/j) , A(Q 2 ,/j) , (7 s 
Q 2 Q A Q6 + -' [) 

where the coefficients Ad>2 contain condensate densities with growing dimensions. In 
the above, \x is the ultraviolet renormalization scale in the loop diagrams. The OPE is 
applicable at sufficiently large Q 2 , provided that the coefficients Ad are proportional to the 
powers of the light-quark masses and/or to the condensate densities, the latter being of 
O(Aqcd) in some power. Hereafter we neglect the u, d- quark masses with respect to m s . 



2 A possibility to assess higher Gegenbauer moments is provided by the model of nonlocal conden- 
sates [Ml. 
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Figure 1: Diagrams contributing to OPE of the correlation function (TJJ) in leading or- 
der: upper row: perturbative loop and quark-condensate diagrams, lower row: examples 
of the gluon-, quark-gluon- and four- quark- condensate diagrams. The wavy (curly) lines 
denote external currents (gluons), the lines with crosses are vacuum fields. The shaded 
oval distinguishes the current with the derivative. 

In the expansion (J7J), the terms proportional to 1/Q d>6 are also neglected, since already 
the contribution of the d = 6 term is quite small. Hence, only the vacuum condensates 
with dimension d < 6 are taken into account. 

The gluon radiative corrections can be systematically included in each term of the 
OPE (J7J). Diagrammatically, the a s - (a 2 -) corrections correspond to all possible one-gluon 
(two-gluon) insertions in the Fig. 1 diagrams. The second small parameter entering the 
OPE of the correlation function ([6]) is the ratio m 2 s /Q 2 . Hence, each of the coefficients 
in (J7J can be cast into a form of a generic double expansion: 

mq\») = 4 0>0) + ) 4 1,0) + 0) + (g 4 M) 

, ( m l\ 2 n (0,2) ,(<*s\ fm 2 s \ (i,i) ( , 

where the coefficients multiplying {a s /^) l {m 2 s /Q 2 ) k depend on ln(yU 2 /Q 2 ) = Iq- 

It is important to assess the numerical role of the small parameters in the combined 
expansion ©. We expect to use OPE at Q 2 ~ 1 GeV 2 . Taking a s (l GeV) = 0.47 [15] and 
a conservative upper limit for the strange quark mass, m s (l GeV) < 150 MeV, one has 
m2 s /Q 2 < 0.02 <C ol s /ix ~ 0.15. Hence, the perturbative 0{a k s ) contributions to the OPE 
are expected to be more important than the 0((m 2 /Q 2 ) k ) terms with the same power. In 
particular, in the first line of (jSJ) the second-order, a 2 - correction is expected to be of order 
of the m 2 /Q 2 -term. Moreover, the observed hierarchy allows one to neglect all "mixed" 



4 



0(a l s (m 2 /Q 2 ) k ) terms with i, k ^ 0. In what follows, we restrict ourselves to the previously 
known first-order corrections in m 2 /Q 2 in the d = 2,4 terms of OPE, neglecting them in 
the d = 6 term. 

So far only the 0(a s ) correction to the quark-condensate contribution A± was calcu- 
lated [31 H], that is, the coefficient in (jSJ). Here we repeat this calculation and, in 
addition, compute the 0(a s ) and 0(a 2 ) terms in A2 and the 0(a 2 ) term in A± that is, the 
coefficients a^ 10 '* , a^ 20 " 1 and af°\ respectively, in (JSj) . Hence, for the largest d — 2,4 terms 
of the OPE ([7]) the NNLO accuracy in a s is achieved. For the subleading d = 6 term in 
the OPE we retain the known LO result. 

The gluon radiative corrections are computed in the standard MS-scheme of renormal- 
ization. The well developed techniques of loop calculations are employed, in particular, 
the programs QGRAF [17J, FORM pj, and MINCER [151120]. The following results are 
given for rif = 3. Including the new 0(az s ) and 0(a 2 ) contributions, we obtain the d = 2 
term in (JTj) originating from the perturbative contribution: 



4vr 2 



1 + 



7T 



26 10, 

1 la 
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7T / 
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14449 , 
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The d — 4 contribution in (JTj) generated by the quark-condensate term has the following 
expression to 0(a 2 ): 



-m s (ss) 1 - 



at. 



TT 



112 8 
~27 + 9 Q 



28135 A ., n . 218, 49, 2 

4C(3) -I In H l n 

1458 SV ; 27 Q 81 Q 



+ 2 



Q 2 



- mjjiii) [ ^ + 



a. 



7T 



59 49, 



(10) 



where (gg) = (0|gg|0), (q = s,u) is the quark-condensate density. 

Finally, the d = 6 term in ([7]) contains the LO quark-gluon, gluon and four-quark 
condensate contributions [H H] : 



2 1 32 

A(<2 2 ,A0 = + -m^G 2 ) (1 + l Q ) - ^a s [ (ss) 2 - {uuf 



[ID 



where (sGs) = (0\sa II1/ g s G a ^{X a /2)s\0), (G 2 ) = (0\^G a ^G aiMU \0) , and the four-quark 
condensate is factorized into the products of two quark condensate densities, assuming 
isospin symmetry (dd) = (uu). 
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Comparing with the previous calculations [31 H], we observe a difference in the 0(a s ) 
s-quark condensate term. More specifically, the constant term 112/27 in the first line of 
( flOl) has to replace 124/27 in the corresponding expressions obtained in [3 J and [4]H This 
difference has, however, only a minor influence on the numerical results. 

The form of the coefficient function multiplying rn^/Q 2 in (J9j) and m 2 s (G 2 ) in (llip 
deserves a separate discussion. One has to emphasize that, in order to achieve a com- 
plete factorization of small and large scales in OPE, all logarithms of the small parameter 
m s have to be removed from the coefficient functions leaving only the powers of m 2 /Q 2 . 
This procedure, understood long ago [ZD, E2J [231 EU EH], generates terms proportional 
to ln(fi 2 /Q 2 ) instead of ln(m 2 /Q 2 ), e.g., in the LO parts of Ai and Aq in (jHD and ( fill , 
respectively. These logarithms were not properly treated in previous calculations. Note 
that a simple replacement \n(m 2 /Q 2 ) — > ln(/i 2 /Q 2 ) in the coefficient functions can miss a 
constant term which has to be added to the logarithm. Let us, for example, explain the 
calculation of the coefficient function in (llip . taking into account the mixing of (sGs) and 
m s (G 2 ) terms under renormalization in MS scheme. 

We isolate the contributions of the quark- ant iquark-gluon and gluon condensates to the 
correlation function ([6]) and write them in the following convenient form: 



n(Q 2 



3Q 6 



C 1 {Q 2 )m s {sGs) + -C 2 (Q 2 )m 2 s (G 2 ) )+... 



(12) 



where all other contributions indicated by ellipses are not important for this discussion. 
The two terms shown in (fl2"j) are generated by the renormalized d = 6 operators m s sGs 
and m 2 G 2 . To calculate the coefficient functions Ci^{Q 2 ) at D = 4 — 2e it is sufficient to 
write the pattern of the mixing of these two operators in the form 



rriesGs 



m s sGs H m 2 G 2 

2e s 



(13) 



where the index "nr" indicates that the operators on r.h.s. are constructed from bare, non- 
renormalized quark and gluon fields. The Z factor of the multiplicative renormalization of 
{sGs) nr can be put to unity in this approximation. To proceed, we need the expansion of 
the operator product in before the vacuum average is taken: 



z / d 4 x e iq x T{u(x) W s(x), s(0)wD a m(0)} 



3Q 6 



m s sGs H m 2 G 2 

2e s 



+ \c 2 (Q 2 ) 



mlO 1 



+ ... 



(14) 



where we use (JT3"]) . The coefficient functions Ci^iQ) can now be obtained by projecting 
both parts of ( fl4l) onto the suitable free quark-gluon states. Sandwiching ( !T4"|) between the 



3 The agreement between [3] and our result is restored when the renormalization procedure for the 
one-loop diagrams used in [3] is corrected by taking into account the 0(e)-term arising in D = 4 — 2e from 
the tree-level diagram with quark condensate. 
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(sg\ and \s) states, and computing the relevant diagrams we find: 

G 1 = l + ' 7 -e. (15) 
o 

Note that one has to retain the 0(e) term for the finite tree- level diagrams. Furthermore, 
projecting ffT4|) onto the two-gluon state, we obtain: 

1 1 13 

-Ca + C 2 = - + — + l Q . (16) 

e e o 

After substituting G\ from ( fT5j) to ( [IB"]) , the 1/e poles cancel each other and the desired 
results for the renormalized coefficient functions are obtained at e — > 0: 

Ci(Q 2 ) = l, C 2 (Q 2 ) = l + l Q . (17) 

We again emphasize that one has to keep the 0(e) term in the expression for the tree-level 
coefficient function C\(Q 2 ) in order to get the non-logarithmic contribution to C 2 (Q 2 )- 
The absence of 1/e poles in the above calculation can also be interpreted as a result of 
the cancellation between ultraviolet and infrared divergences. Indeed, the 1/e pole in the 
mixing of operators emerges as an ultraviolet divergence, whereas 1/e in the loop diagrams 
used to calculate the coefficient functions has an infrared origin. A similar derivation is 
applied to the mixing of 0(m s (ss)) and 0(m^) terms in OPE, yielding the constant term 
5/2 that accompanies hi(p, 2 /Q 2 ) in the coefficient function (Q. 



4. After the correlation function is calculated, the derivation of the sum rule follows 
the standard procedure. The OPE (j7|) with the coefficients given in f[9l. f[TUj) and ffTTj) is 
equated to the hadronic dispersion relation 

aiq ^Mit 2 + J dB m. (18) 

m K ~ Q J s ~ 1 

Sh 

In the above, possible subtractions are ignored in anticipation of the Borel transformation. 
The residue of the kaon pole is obtained by combining the matrix element (j3J) with the 
definition of the kaon decay constant. 

The spectral density p h (s) includes the contributions of hadronic continuum and reso- 
nances with J p = 0",1 + and strangeness: Ktttt, K*tt, Kp, Ki(1270), i^i(1400),... . Ac- 
cordingly, the lower limit of integration is Sh = (ttlk + 2m,,-) 2 , the invariant mass squared 
of the lightest continuum state in this channel. To approximate p h (s), we employ the 
quark-hadron duality approximation: 

P h (s)e(s-s h ) = P OPE (s)e(s-s«), (19) 
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where sjf is the effective threshold, and the spectral density: 



„ope. 



s, /i) = — ImII(s, fi) 

7T 

| m s (ss) f 8a, 



10a,. 



+ 



9tt 



7T 



AtT 2 S \ 97T 

2 r218 98, 

1 I 

L 27 81 *J 



7T 



2 r 14449 605, 
+ —L 



972 



162 J 



3m 2 



m s (uu) 49 fa; s \ 2 m 2 (G 2 ) 
ix ) 3s 3 



(20) 



with / s = ln(// 2 /s), is obtained by calculating the imaginary part of H(Q 2 ) at positive 
s = —Q 2 . In the above, the running parameters (a s , m s ) are taken in the MS scheme and 
normalized at the scale /i. 

The next step is the Borel transformation of (|T8l) . which eliminates the subtraction 
terms and suppresses the integral over p h (s) , so that the resulting relation becomes less 
sensitive to the duality approximation. The transformed invariant amplitude IT (g 2 ) has 
the following form: 



n(M 2 



+ 



+ 



4tt 2 
605 



1 + — 

7T 



26 10,, 



((/M + 7,) 2 -y)-fc(3) 



324 

m s {ss) 
M 2 



7T 



366659 14449 



11664 



972 



(Im + 1e) 



M 2 



2+Im + 1e 



I _ ^1 



27 + g(^M + 7e) 



7T 



18127 556,, . 



m s {uu) ( Aa s | 
M 2 



4C(3) 



+ 



M 2 



9tT \ 7T 



79 49 „ 

162 + 81 (/M + 7s) 



m s (sGs) + m 2 (G 2 ) 



3M 4 



6M 4 



~ + Im + 1e 



167ra s 
27M 4 



(ss) 2 - {uu} 2 



(21) 



where Im = ln(/i 2 /M 2 ) and 7^ is the Euler constant. Finally, the sum rule for the kaon 
Gegenbauer moment af obtained from ( fT8l) using ( [19]) reads: 



.A' 



3/1 



(00 
n(M 2 ) - y rf sp OPS ( 



s e 



s/M 2 



(22) 



where the functions I1(M 2 ) and p (s) are given in (|2ip and fl20l) . respectively. An equiv- 
alent form of the sum rule where the entire r.h.s. is represented as a duality integral over 
the spectral density, was used in [SIS]. 
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5. To perform the numerical analysis of the sum rule (1221) . we specify the relevant 
input parameters, starting from the kaon mass m^- = 493.58 MeV and decay constant 
fx = 159.8 ± 1.4 ± 0.44 MeV [16] . For the strange quark mass we adopt: 

m s (2 GeV) = 98 ± 16 MeV (23) 

which covers the intervals m s (2 GeV) = 105 ± 6 ± 7 MeV [9] (with uncertainties added in 
quadrature) and m s (2 GeV) = 92 ± 9 MeV [10] from the recent QCD sum rule determina- 
tions with 0(a 4 ) accuracy. 

The running of the coupling a s (mz) = 0.1176 ± 0.002 [16] to the lower scale /i is taken 
with the 4-loop accuracy applying the program presented in [26] . The default value of the 
renormalization scale is /i = 1 GeV, where the expansion parameter is a s (l GeV)/7r = 
0.15 ± 0.01 and m s (l GeV) = 128 ± 21 MeV is obtained from (123]). 

The nonstrange quark condensate density (qq) (q = u, d) is fixed from the Gell-Mann- 
Oakes-Renner relation (qq) = —m 2 K f 2 / ] (2(m u +m ( i)\. Taking the non-lattice averages for the 
u- and <i-quark masses, m u (2 GeV) = 3±1 MeV and md{2 GeV) = 6.0 ±1.5 MeV, together 
with f„ = 130.7 ± 0.1 ± 0.36 MeV and = 139.57 MeV from [16], we get (qq(2 GeV)) = 
-(0.264+°;°^ GeV) 3 , and correspondingly, (qq(l GeV)) = -(0.242±°;° 2 * GeV) 3 . As men- 
tioned above, in all other formulae the light-quark masses are neglected. Furthermore, the 
ratio of strange and nonstrange condensates (ss) / (qq) = 0.8 ± 0.3 is adopted [271 ES 129] . 

The accuracy of the condensate densities with higher dimensions is less important. In 
particular, we take for the quark-gluon condensate density the standard parameterization 
(sGs) = ml(ss)(l GeV) with m 2 = 0.8 ± 0.2 GeV 2 [301 129] (neglecting the running), and 
a very wide interval for the gluon condensate density (G 2 ) = 0.012 ± 0.012 GeV 4 . We also 
allow for a varying factor 0.1 -r 1 multiplying the square of the quark condensate density 
in the factorization relation for the four-quark condensates. 

The Borel parameter interval M 2 = 1.0 — 2.0 GeV 2 is adopted, so that the d = 6 
contribution in I1(M 2 ) remains lower than 30%. The upper limit is taken to avoid too large 
contributions of the excited and continuum states estimated with duality approximation. 
The dependence of af- on M 2 for Sq = 1.05 GeV 2 and at central values of all input 
parameters is plotted in Fig. 2, showing a remarkable stability. There we also plot separate 
contributions to the sum rule originating from the d — 2,4, 6 terms of the expansion ([7]). 
The threshold parameter Sq = 1.05 GeV 2 was fixed in [3] from the sum rule for f^. 
In fact, the -dependence of the sum rule result turns out to be rather weak, if one 
varies the threshold in rather wide limits Sq = 0.9 — 1.4 GeV 2 , the dependence of 
on M 2 remaining flat. The reason is that the spectral density f)20p is small in the region 
above the duality interval for the kaon. Therefore, although the pattern of hadronic states 
(resonances and continuum states) in this region is rather complicated, we expect that the 
integral over these states is also small. Hence the sum rule result will not noticeably change 
if one modifies the duality ansatz for the hadronic spectral density, e.g. by adding more 
resonances and correspondingly increasing the threshold parameter (as shown e.g., in [3]). 

The numerical prediction of the sum rule is represented in the form: 

af (1 GeV) = 0.100 ± 0.003| S r ± 0.003| Qs ± 0.035| ms ± 0.022| mq ± 0.013| cond , (24) 
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Figure 2: The kaon Gegenbauer moment af (1 GeV) calculated from the QCD sum rule as 
a function of the Borel parameter (solid); the contributions of d = 2, d = 4 and d = 6 
terms to |7|) are shown with dashed, dotted and dash-dotted lines, respectively. 



where the first ( "SR" ) error is the combined uncertainty of of due to the variation of M 2 
and Sq , representing a sort of an "intrinsic " uncertainty of the sum rule. The subsequent 
errors correspond to the individual variations of a s (mz), m s (2 GeV) and m u ^{2 GeV) 
(that is, (qq)), within the adopted intervals. The last error ("cond") shown in (T24"j) is a 
combined uncertainty due to variation of all remaining condensate parameters. Finally, 
adding the individual uncertainties in quadrature we obtain the interval ([T]). 

Let us now discuss the structure of the perturbative series as it follows from the numer- 
ical analysis. First, we observe that gluon radiative corrections substantially enhance the 
perturbative d = 2, 0(m 2 ) term and suppress the quark condensate term. Numerically, 
for the 0{m 2 ) term in (I2TI) at // = M = 1 GeV the following a s -expansion is obtained (in 
MS-scheme): 



n K)(i GeV) 



1 GeV) 



4vr 2 



3.53 



a s {\ GeV) 



7T 



+ 33.7 



a s (l GeV) 



7T 



(25) 



revealing a poor convergence in a s at a s (l GeV) /n ~ 0.15. It is not surprising, if one 
recalls the higher-order perturbative corrections calculated for the m s determination from 
r decays or from QCD sum rules based on scalar /pseudoscalar correlation functions, where 
the situation is in fact similar. For instance, in the 0{m 2 s ) part of the Borel-transformed 
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scalar /pseudoscalar correlation function [8j [9] one has (see e.g., eq.(16) in [9]): 



, a s (l GeV)\ faJl GeV)\ 2 
1 L..S21 ( — J +21.98 — -j 



7T 



+ 53.1| a ' (1 ° eV) ) 3 + 31.6( a - ( \ GeV) ) 



3m^(l GeV) 



1 + 0.72 + 0.49 + 0.18 + 0.02) 



(26) 



The ratio of (I26I) (taken with O(a^) accuracy) and fT25|) has a much better convergence: 

2 



n (5)"(<)(i GeV) _ 3 
lK m i)(l GeV) ~ 2 



i + 1 .29i"- (1 ^ v) )- 1 6.2("-( 1 f v )) + 



(27) 



where the corrections cancel each other to a great extent. 

In the quark condensate contribution, the radiative corrections are less sizable: at 
fi = 1 GeV the numerical hierarchy of the terms multiplying the strange quark condensate 
density in ( fTUI) reads: 



n (-(-))(i GeV) = m s (ss) ( 1 - 3.77 ('^L^XT) _ i . 8 



(2* 



where we again find that the second-order correction in a s is important. 

We conclude that after taking into account the NNLO perturbative corrections, the 
numerical pattern of the sum rule for drastically changes: the coefficient of the d = 2 
term gets enhanced, whereas the d = 4 term decreases. Note also that including higher- 
order perturbative corrections in the sum rule, makes more consistent the use of the input 
parameters, such as m s , determined with a high accuracy, up to 0(af). 

Finally, the d = 6 subleading contributions to the sum rule play an important role 
in providing the Borel stability. The mixed quark-antiquark-gluon condensate dominates 
numerically in (TTTT) yielding a negative contribution to the Borel-transformed correlation 
function and stabilizing the whole sum of OPE terms. Computation of the radiative correc- 
tions to this term is a difficult task beyond our scope. Having in mind the large uncertainty 
of the mixed condensate density, we expect that the a s -corrections to the coefficient func- 
tion of the (sGs) term will hardly improve the overall accuracy of its contribution. 

Finally, we checked that due to the enhanced precision in a s , the dependence of the sum 
rule prediction for af on the renormalization scale /i becomes small. Taking fi > 1 GeV 
we calculate af (1 GeV) by rescaling the sum rule result with the NLO scale dependence 
([5]). The LO and NLO logarithmic dependences naturally cancel out in d = 2,4 terms, 
leaving a very mild residual scale-dependence due to the unaccounted a s -correction to the 
(sGs) term. 
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6. Concluding, we have calculated the NNLO gluon radiative corrections to the QCD 
sum rule for the first Gegenbauer moment of the kaon distribution amplitude. The correc- 
tions turned out to be numerically important, they change the relative magnitude of the 
d = 2 (loop diagrams) and d — 4, 6 (condensate) terms in the OPE, improving also the 
Borel stability of the sum rule. 

The uncertainty of af is still large and amounts up to 40%, due mainly to the limited 
precision of the light quark masses: m s directly entering the sum rule and m u ^ determining 
the quark-condensate densities via Gell-Mann-Oakes-Renner relation. A better determi- 
nation of the ratio of strange and nonstrange condensates and of the mixed quark-gluon 
condensate are another possibilities of reducing the theory error. The calculation of the 
radiative corrections to the quark-gluon condensate contribution can also improve the ac- 
curacy. This however requires a dedicated computational effort including an analysis of the 
whole basis of the dimension-six operators. The weak dependence of the sum rule on the 
threshold parameter Sq indicates that the one-resonance (kaon) duality ansatz is quite sat- 
isfactory. Still, an additional analysis, including the axial-vector and radially excited kaon 
resonances could provide a better understanding of the duality pattern in this channel. 

Our result for af" is somewhat larger than the previous estimates [31 HI [7], and the 
uncertainty, we believe, is more realistic. For comparison we also quote the two recent 
lattice QCD determinations of this parameter: af (2 GeV) = 0.0453 ± 0.0009 ± 0.0029 [31] 
and af(2 GeV) = 0.048 ± 0.003 |32j, which have achieved a rather small error. By 
evolving our result (CQ) to this scale with the help of the scale- dependence ([5]) we find 
af(2 GeV) = 0.08 ±0.04, that is, within uncertainties, only a marginal agreement with the 
lattice results. Let us note that both lattice determinations use a linear extrapolation in 
m s (kaon mass squared) inspired by ChPT in the leading order (see also [33]). We would 
like to stress again that in our analysis the contribution to the sum rule proportional to 
w? s is enhanced while the term proportional to m s (ss) is suppressed, both enhancement 
and suppression being caused by the radiative corrections. In the language of ChPT, this 
observation could indicate an important role of the next-to-leading terms in the expansion 
of di in the kaon mass. 
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